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The  mechanical  response  of  p-phenylene  terephthalamide  (PPTA)  single  fibers  when  subjected  to  uniaxial  compression  is 
investigated  computationally  using  coarse-grained  molecular  statics/dynamics  methods.  In  order  to  construct  the  coarse-grained 
PPTA  model  (specifically,  in  order  to  define  the  nature  of  the  coarse-grained  particles/beads  and  to  parameterize  various 
components  of  the  bead/bead  force-field  functions),  the  results  of  an  all-atom  molecular-level  computational  investigation  are  used. 
In  addition,  the  microstructure/topology  of  the  fiber  core,  consisting  of  a  number  of  coaxial  crystalline  fibrils,  is  taken  into  account. 
Also,  following  our  prior  work,  various  PPTA  crystallographic/topological  defects  are  introduced  into  the  model  (at  concentrations 
consistent  with  the  prototypical  PPTA  synthesis/processing  conditions).  The  analysis  carried  out  clearly  revealed  (a)  formation  of 
the  kink  bands  during  axial  compression;  (b)  the  role  of  defects  in  promoting  the  formation  of  kink  bands;  (c)  the  stimulating 
effects  of  some  defects  on  the  fiber- fibrillation  process;  and  (d)  the  detrimental  effect  of  the  prior  compression,  associated  with 
fiber  fibrillation,  on  the  residual  longitudinal-tensile  strength  of  the  PPTA  fibers. 


1.  Introduction 

Fibers  of  p-phenylene  terephthalamide  (PPTA),  available 
commercially  as  Kevlar,  Twaron,  Technora,  and  so  forth,  are 
characterized  by  high  specific  axial  stiffness  and  high  specific 
tensile  strength.  These  fibers  are  often  referred  to  as  ''ballistic 
fibers''  since  they  are  commonly  used  in  different  ballistic- 
and/or  blast-protection  systems.  Such  systems  are  required  to 
provide  a  high  level  of  penetration-resistance  against  large- 
kinetic- energy  projectiles,  such  as  bullets,  detonated-mine- 
induced  soil  ejecta,  fragments/shrapnel  from  improvised 
explosive  devices  (lEDs),  and  shells  or  turbines.  To  satisfy 
this  requirement,  PPTA  fibers  are  normally  used  as  (a) 
filament  constituents  in  twodimensional  or  threedimensional 
woven-fabric  protective  structures  (e.g.,  "bulletproof  vests"); 
(b)  reinforcements  in  high-performance  (typically,  polymer- 
matrix,  "rigid- armor")  composites;  or  (c)  energy-absorbing 
reinforcements  in  resin-starved  flexible  armor-grade  com¬ 
posites  [1-4]. 


Given  the  complex  behavior/functionality  of  the  afore¬ 
mentioned  protection  systems,  their  design  and  development 
increasingly  involve  the  use  of  computer-aided  engineering 
(CAE)  methods  and  tools.  In  turn,  the  latter  employ  material 
models  to  describe  the  behavior  of  fibers  and  structures  under 
high-rate  loading  conditions.  With  the  utility  of  the  CAE 
methods  and  tools  largely  determined  by  the  fidelity  of  these 
models,  their  development  must  account  for  the  hierarchi¬ 
cal/multiscale  architecture  of  the  fibers-based  structures,  as 
well  as  of  the  fibers  themselves  [5-8]. 

The  present  work  deals  with  the  compressive  axial  behav¬ 
ior  of  single  PPTA  fibers  and  the  effect  of  this  behavior 
on  the  global/macroscopic  response  of  the  ballistic-  and 
blast-protection  systems.  Consequently,  the  main  topics 
overviewed  in  the  remainder  of  this  section  include  (a)  the 
basics  of  PPTA  fibers;  (b)  PPTA  microstructural  hierarchy; 
and  (c)  overview  of  the  recent  materialmodels  for  PPTA  fiber- 
reinforced  polymer-matrix  composites. 
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1.1.  The  Basics  ofPPTA  Fibers 

1.1.1.  Classification.  PPTA  fibers  can  be  categorized  as 
polyamide-type  (polymeric)  fibers.  Polyamides  are  generally 
divided  into  two  main  classes:  (i)  nonaromatic  polyamides 
(e.g.,  nylon-6,6)  and  (ii)  aromatic  polyamides  or  aramids 
(e.g.,  Kevlar,  Twaron,  Technora,  etc.).  It  is  the  aromatic 
polyamides  which  are  the  subject  of  the  present  work. 

1.1.2.  PPTA  Molecules/ Chains.  The  PPTA  repeat  units/mers 
appear  in  two  distinct  (a)  trans  Con  the  opposite  side''); 
and  (b)  cis  Con  the  same  side")  molecular-level  confor¬ 
mations.  These  conformations  are  shown  using  a  ball-and- 
stick  representation  in  Figures  1(a)  and  1(b),  respectively. 
To  help  with  interpretation  of  the  molecular  structure,  the 
chemical  elements  are  labeled  in  Figure  1(a).  Examination  of 
Figures  1(a)  and  1(b)  reveals  that  the  repeat  unit  consists  of 
two  phenylene  rings/moieties  (C6H4RR^)  connected  by  an 
amide  linkage  (-CO-NH-).  Although  the  PPTA  molecule 
can  assume  either  the  trans  or  cis  stereo-isomeric  con¬ 
formations,  the  former  conformation  is  almost  exclusively 
observed.  The  reason  for  this  predominance  is  that  the  trans 
conformation  promotes  the  stretching/extension  of  PPTA 
molecules,  lowering  the  system  energy.  In  addition,  extended 
molecules  can  readily  align,  forming  sheets  consisting  of 
nearlyparallel  regularlyspaced  PPTA  chains.  Formation  of 
these  sheets  is  associated  with  further  reductions  in  the 
system  energy.  Regular  stacking  of  these  sheets  leads  to  the 
formation  of  low-energy  crystalline  PPTA  fibrils  (subfiber 
units,  elongated  and  coaxial  with  the  fibers).  In  sharp  con¬ 
trast,  the  cis  conformation  hampers  the  formation  of  these 
low-energy  structures,  which  explains  its  rare  occurrence. 

1.1.3.  PPTA  Sheets.  Due  to  a  large  difference  in  electronegativ¬ 
ity  between  oxygen  and  hydrogen,  the  electron  charge  density 
within  the  PPTA  unit  cell  is  increased  around  the  oxygen 
atoms  and  lowered  around  the  hydrogen  atoms,  giving  rise  to 
the  formation  of  (relatively  strong)  dipoles.  The  presence  of 
dipoles  along  with  (i)  the  lone  electron  pairs  on  the  oxygen 
atoms;  (ii)  close  proximity  of  the  oxygen  and  hydrogen 
atoms  of  neighboring  parallel  PPTA  molecules  within  the 
same  sheet;  and  (iii)  the  operation  of  London-dispersion 
effects  typically  results  in  the  formation  of  (strong)  hydrogen 
bonds.  Within  these  bonds,  the  N-H  species  acts  as  a  bond 
donor  while  O  acts  as  a  bond  acceptor.  The  formation  of 
the  aforementioned  sheet-like  structures  is  believed  to  be  the 
result  of  strong  lateral  hydrogen  bonding  between  adjacent 
parallel  PPTA  molecules. 

1.1.4.  PPTA  Fibrils.  As  mentioned  previously,  stacking  of 
sheets  results  in  the  formation  of  crystalline  fibrils,  the  axes 
of  which  are  aligned  with  the  fiber  axis.  However,  orientation 
of  the  fibril  crystal  structure  within  a  plane  normal  to  the  fiber 
axis  varies  among  the  fibrils.  In  other  words,  PPTA  fibrils 
within  a  single  fiber  resemble  columnar  grains  observed 
in  directionally  solidified  as-cast  metallic  microstructures. 
However,  in  contrast  to  the  columnar  metallic  grains,  the 
PPTA  fibrils  are  generally  found  to  possess  a  superstructure 


consisting  of  small- amplitude  (ca.  20  nm),  high  wavelength 
(ca.  300-500  nm)  accordion -style  "fileats."  The  presence  of 
these  pleats  is  often  cited  as  one  of  the  causes  for  the  supe¬ 
rior  tensile  ductility  of  PPTA  fibrils/fibers  (e.g.,  [9]).  While 
intrasheet  bonding  is  of  a  strong  hydrogen -bonding  type, 
intersheet  bonding  involves  weak  van  der  Waals  interactions 
and  the  formation  of  Ti-type  molecular  orbitals.  Thus,  PPTA 
fibrils  possess  a  substantial  level  of  transverse  anisotropy. 

As  mentioned  previously,  intersheet  bonding  within 
PPTA  fibrils  is  relatively  weak,  and  consequently,  fibrils  tend 
to  possess  various  crystallographic  and  topological  intersheet 
defects  (e.g.,  local  delaminations  [10]).  This  is  the  reason 
that  PPTA  fibrils  tend  to  buckle  during  axial  and  lateral 
(sheet-edge-on)  compressive  loading  and  to  possess  relatively 
inferior  axial/transverse  compressive  properties. 

1.1.5.  PPTA  Fibers.  As  mentioned  previously,  PPTA  fibers 
are  essentially  assemblies  of  parallel  fibrils  (nearly  coaxial 
with  the  fiber  itself).  While  it  is  generally  believed  that  fibril 
bundles  within  the  fiber  are  held  together  by  weak  van  der 
Waals  or  Coulomb  forces,  it  is  possible  that  chain  crosslinking 
and  the  presence  of  tie  fibrils  can  play  a  role  in  the  formation 
of  fibril  bundles.  A  schematic  of  the  basic  architecture  of  the 
single  PPTA  fiber  and  its  relationship  to  the  constituent  fibrils 
is  depicted  in  Figure  2.  Examination  of  this  figure  reveals 
that  the  fibril  bundles  form  the  fiber  core  and  that  the  fibers 
also  contain  a  relatively  thin  skin  region.  While  this  region  is 
comprised  of  highly  oriented  PPTA  chains,  it  is  usually  of  a 
noncrystalline  character. 

1.2.  PPTA  Microstructural  Hierarchy 

1.2.1.  Structural/ Behavioral  Complexity  of  Protective  Systems. 
The  work  presented  in  this  paper  is  directly  concerned 
with  PPTA  single  fibers,  or  more  specifically,  with  the 
axial  compressive  response  of  these  fibers.  However,  the 
broader  goal  of  the  present  work  is  to  help  improve  fidelity 
and  accuracy  of  the  computational  models  for  PPTA-fiber- 
reinforced  polymer-metal  composites  used  in  the  aforemen¬ 
tioned  blast-Zballistic-protection  systems.  These  protective 
systems  generally  possess  complex  morphology  and  behavior 
[5-8],  which  is  mainly  attributed  to  (a)  their  hierarchical, 
multilength-scale  architecture;  (b)  their  nonlinear  and  path- 
and  rate-/time-dependent  mechanical  response;  and  (c)  the 
operation  of  various  phenomena/processes  (e.g.,  fibril/fiber 
damage,  highly  localized  deformation  and  the  associated 
formation  of  kink/deformation  bands,  filament  twisting, 
interfilament  friction  and  sliding,  etc.). 

1.2.2.  Major  Length  Scales.  By  carrying  out  a  detailed  analysis 
of  the  structural  hierarchy  of  the  PPTA-fiber-reinforced 
polymer-matrix  composite  protective  structures,  eight  dis¬ 
tinct  length  scales  were  identified  in  our  recent  work  [5,  6]. 
The  microstructures/architectures  of  each  of  these  length 
scales  are  schematically  depicted  in  Figures  3(a)-3(h).  In  each 
case,  labels  are  used  to  identify  the  main  microstructural 
constituents.  Since  a  detailed  description  of  these  length 
scales  can  be  found  in  [5,  6],  they  will  not  be  repeated  here. 
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(a)  (b) 

Figure  1:  (a)  Trans-  and  (b)  c/s-molecular  conformations  in  typical  PPTA-based  polymeric-material  chains/molecules. 


Figure  2:  Topological  relationship  between  a  single  fiber  and  its 
constituent  fibrils. 


1.3.  Material  Models  for  PPTA  Fiber- Reinforced 
Polymer-Matrix  Composites 

1.3.1.  CAE-Based  Design/Development  of  Protective  Systems. 
Ballistic/blast  protection  systems  based  on  PPTA  fibers  have 
traditionally  been  designed  and  developed  using  legacy 
knowledge,  empiricism,  and  extensive  'fabricate- and-tesC 
procedures.  Given  the  higher  costs  and  often  significantly 
longer  lead  times  associated  with  this  approach,  it  has 
gradually  been  complemented  (and  partly  replaced)  by  the 
appropriate  cost-  and  time-efficient  CAE  analyses.  Recent 
developments  in  the  numerical  modeling  of  transient  non¬ 
linear  dynamics  phenomena,  such  as  those  accompanying 
blast-  and  ballistic-loading  conditions  (e.g.,  [11,  12])  have 
further  accelerated  the  use  of  CAE  analyses  in  the  design  and 
development  of  protective  structures.  However,  the  accuracy 
of  these  CAE-based  results  depends  on  the  fidelity  of  the 
constitutive  models  for  the  PPTA-based  materials  used  in 
constructing  protective  structures  (e.g.,  [13]).  As  demon¬ 
strated  in  our  recent  work  [7],  the  PPTA  material  models 
generally  used  suffer  from  a  number  of  deficiencies.  These 
deficiencies  are  associated  with  (a)  inability  of  the  models  to 


fully  account  for  the  material  response  to  high-deformation- 
rate,  large-strain,  high-pressure  loading  conditions  (typically 
encountered  during  projectile-  or  blast-impact  events);  (b) 
disregard  for  the  effects  of  prior  treatment/handling  of  the 
constituent  fiber  materials  on  the  resulting  material  prop¬ 
erties;  and  (c)  a  lack  of  inclusion  of  the  changes  in  the 
material  properties  brought  about  by  the  protective-structure 
manufacturing  process  chain. 

1.3.2.  Recent  Improvements  in  PPTA  Material  Model.  In  the 
same  work  [7],  it  was  demonstrated  that  the  first  step  towards 
removing  these  deficiencies  is  to  identify  and  characterize  key 
phenomena  and  processes,  occurring  at  different  length/time 
scales  and  to  quantify  their  contributions  to  the  overall 
behavior/properties  of  the  subject  (continuum  type)  material. 
In  this  regard,  the  PPTA  material- constitutive  models  were 
improved  in  our  recent  work  [5-8]  as  follows:  (a)  inclu¬ 
sion  of  the  effect  of  various  fiber-  and  subfiber-length-scale 
microstructural  features  (including  the  defect  structure)  on 
the  material  dynamic  response  [5,  6];  (b)  incorporation  of 
the  inelastic  deformation  effects  associated  with  transverse 
loading  and  their  coupling  to  the  longitudinal  mechanical 
properties  of  the  PPTA  fibers  [8];  and  (c)  accounting  for 
the  effects  of  prior  longitudinal-compressive  and/or  axial- 
torsional  loading  on  the  tensile  behavior  of  the  PPTA  fibers 
and  so  accounting  for  the  effect  of  fiber  handling  during 
PPTA-structure  fabrication  [14]. 

1.3.3.  Proposed  Additional  Improvements.  In  the  present 
work,  an  attempt  is  made  to  further  improve  the  PPTA- 
material  model  used  in  the  CAE  analyses.  Specific  improve¬ 
ments  deal  with  the  contribution  of  the  axial  compressive 
response  of  single  PPTA  fibers,  which  results  in  the  formation 
and  propagation  of  kink  bands  (fiber-surface  irregularities 
in  the  form  of  disclinations/protrusions),  Eigure  4.  Since 
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Figure  3:  Microstructural  hierarchy  consisting  of  eight  distinct  length  scales  encountered  in  fiber- reinforced  polymer-matrix  composites. 


PPTA  fibers  can  experience  axial-compression  (including  the 
one  accompanying  fiber  bending)  during  fabrication  of  the 
reinforcing  substructures  (e.g.,  during  textile  manufactur¬ 
ing  of  ballistic  fabric),  the  PPTA  fiber- reinforced  polymer- 
matrix  composites  and  the  associated  protective  structures 
typically  contain  kinked  fibers  (fibers  with  inferior  mechan¬ 
ical  properties).  In  addition,  snap-back  of  elastic  unloading 
following  fracture  of  the  stretched  fibers  (the  process  which 
often  accompanies  ballistic-Zblast-impact  loading)  is  often 
accompanied  by  the  formation  of  a  series  of  kink  bands. 
The  resulting  inferior  mechanical  properties  of  the  kinked 
fibers  can  seriously  compromise  multihit  performance  of  the 
protective  structures. 

1.4.  Main  Objective.  The  main  objective  of  the  present  study 
is  to  extend  our  prior  molecular- type  simulation  work  [5- 
8,  14]  in  order  to  (i)  provide  insight  into  the  phenomena 


Figure  4:  A  schematic  of  a  series  of  kink  bands  formed  over  a  short 
segment  of  a  PPTA  fiber. 

and  processes  associated  with  the  formation  and  propagation 
of  kink  bands  within  single  PPTA  fibers  during  axial  com¬ 
pression.  Presently,  it  is  not  clear  if  kink  bands  are  merely 
the  result  of  highly  localized  crystallographic  slip  within  the 
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crystalline/paracrystalline  fibrils  or  are  the  product  of  flow- 
induced  chain/fibril  buckling/localized  bending  (generally 
referred  to  as  ''fibrillation')  [15];  and  (ii)  identify  the  ways  and 
means  for  passing  the  information  regarding  the  operation 
of  (discrete)  fiber  length-scale  phenomena  and  processes 
(e.g.,  formation  and  propagation  of  kink  bands)  to  the 
(continuum)  lamina/laminate  length  scales. 

It  should  be  noted  that  while  the  majority  of  the  bal¬ 
listic/blast  protective  systems  fabricated  today  are  based  on 
polymer-matrix  composites  containing  Kevlar  KM2  fiber- 
based  reinforcing  substructures  (e.g.,  fabric  mats),  the  subject 
of  the  present  work  are  generic  PPTA  fibers.  This  should 
not  be  considered  as  a  major  simplification  since  the  behav¬ 
ior  and  properties  of  various  Kevlar  and  Twaron  grades 
of  fibers  are  generally  found  to  be  dominated  by  their 
main  PPTA  constituent.  Hence,  the  inclusion  of  fiber-grade- 
specific  chemical  and  molecular-microstructure  details  was 
judged  unnecessary  (in  particular,  considering  the  associated 
additional  computational  cost). 

2.  Effect  of  Prior  Compression  on 
Fiber  Tensile  Strength 

As  mentioned  earlier,  kink  bands  resulting  from  the  appli¬ 
cation  of  a  prior  axial/longitudinal  compressive  loading 
can  negatively  affect  axial-tensile  strength  of  PPTA  fibers. 
Examination  of  the  open- domain  literature  revealed  two 
investigations  which  are  particularly  pertinent  to  the  present 
work.  A  brief  overview  of  these  prior  investigations  is  pro¬ 
vided  in  the  remainder  of  this  section. 

2.1.  Dobb  et  al,  1981  [10].  Dobb  et  al.  [10]  carried  out 
longitudinal  compression  tests  on  single  PPTA  fibers  and 
examined  the  fiber  microstructure  before  and  after  testing 
in  order  to  identify  microstructural  changes  responsible  for 
the  inferior  axial-compressive  strength  of  the  fibers.  The 
main  findings  of  Dobb  et  al.  can  be  summarized  as  follows: 
(a)  the  inferior  axial  compressive  strength  of  single  PPTA 
fibers  is  a  mere  consequence  of  the  formation  of  kink  bands 
and  subsequent  fiber  buckling;  (b)  the  aforementioned  fibril 
pleats  appear  not  to  play  a  major  role  in  the  onset  of  kink 
band  formation  or  in  band  propagation;  (c)  formation  of 
kink  bands  has  been  observed  to  be  closely  linked  with 
relatively  weak  van  der  Waals-type  intersheet  bonding,  and 
the  associated  weak  lateral  support  experienced  by  the 
PPTA  sheets.  These  conditions  can  lead  to  the  phenom¬ 
ena  such  as  localized  intersheet  shear,  sheet-delamination 
and  sheet/chain  buckling,  all  of  which  appear  to  stimulate 
the  formation  of  kink  bands;  and  (d)  the  presence  of  the 
crystalline  structure  within  the  PPTA  fibrils  promotes  the 
formation  of  the  kink  bands  through  the  coordinated  motion 
(accompanying  highly  localized  crystallographic  slip)  of  the 
well- organized  PPTA  chain  segments. 

2.2.  Deteresa  et  al.,  1984  [16].  Deteresa  et  al.  [16]  carried 
out  a  systematic  experimental  investigation  of  the  effect  of 
prior  axial  compression  on  the  longitudinal  tensile  strength 


of  single  PPTA  fibers.  The  mechanical  tests  were  supple¬ 
mented  by  optical  and  scanning  electron  microscopy  investi¬ 
gations  in  order  to  relate  changes  in  the  fiber  microstructure 
induced  by  prior  (longitudinal-compressive)  deformation  to 
the  observed  losses  in  the  fiber  longitudinal-tensile  strength. 
The  results  obtained  by  DeTeresa  et  al.  can  be  summarized 
as  follows:  (a)  prior  longitudinal  compression  produces 
regularly  spaced  kink  bands  in  the  single  PPTA  fibers;  (b) 
the  presence  of  these  bands  within  a  single  fiber  reduces 
its  longitudinal- tensile  strength.  Typically,  a  loss  of  ca.  10% 
in  the  fiber  longitudinal-tensile  strength  is  associated  with 
the  prior  application  of  3%  axial- compressive  strain;  and 
(c)  axial  compression  tends  to  induce  well-defined  discrete 
kink  bands,  suggesting  the  dominant  role  of  crystallographic 
slip.  This  behavior  is  in  sharp  contrast  to  the  effect  of  prior 
torsional  deformation,  which  tends  to  induce  diffuse  shear 
bands  and  fibrillation. 


3.  Coarse-Grained  Fiber-Scale 
Computational  Analysis 

As  mentioned  earlier,  coarse-grained  fiber-scale  computa¬ 
tional  methods  have  been  employed  in  the  present  work 
in  order  to  (a)  provide  insight  into  the  phenomena  and 
processes  associated  with  the  formation  and  propagation  of 
kink  bands  within  single  PPTA  fibers  during  axial  compres¬ 
sion  and  (b)  assess  the  role  of  crystallographic/topological 
defects  in  the  kink  band  formation  process  and  in  the  relative 
importance  of  localized  crystallographic  slip  and  fibrillation 
processes. 

Within  the  coarse-grained  computational  framework, 
atoms,  ions  and  interatomic  bonds  are  not  modeled  explic¬ 
itly,  but  rather,  individual  atoms/ions  are  combined  into 
larger  single  particles,  referred  to  as  "beads','  where  each 
bead  represents  the  collective  degrees  of  freedom  of  the 
constituent  atoms/ions.  Neighboring  beads  residing  within 
the  same  chain  are  then  connected  using  bead/bead  bonds 
(also  referred  to  as  connectors),  while  nonbonded  beads  are 
allowed  to  interact  via  van  der  Waals-like  and  Coulomb-like 
force  fields  (the  term  "force  field"  will  be  defined  later).  Since 
the  coarse-grained  computational  framework  employs  larger 
particles  (relative  to  the  particle  sizes  encountered  in  the 
so-called  all-atom  computational  analyses),  larger  computa¬ 
tional  domains  can  be  analyzed  (simulating  the  maximum 
number  of  particles  permitted  by  the  available  computational 
and  storage  resources).  In  addition,  by  integrating  out  the 
high-frequency  vibrations  of  the  constituent  atoms/ions  and 
allowing  the  system  dynamics  to  be  controlled  by  the  com¬ 
pliant/soft,  bead/bead  interactions  (rather  than  by  stiff/hard, 
atom/atom  interactions),  a  significantly  longer  time  scale  is 
attained  (relative  to  that  which  could  be  afforded  by  the  all¬ 
atom  computational  method). 

It  should  be  noted  that  while  only  coarse-grained  molecu¬ 
lar  methods  and  tools  are  used  in  the  present  work,  the  results 
of  the  all- atom  molecular  dynamics  simulations  carried  out  in 
our  recent  work  [5-8, 14]  were  utilized.  Specifically,  all-atom 
molecular  dynamics  calculations  described  in  [5-8, 14]  were 
used  to  (a)  identify  the  nature  and  the  number  of  types  of 
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coarse-grained  beads;  (b)  for  each  bead  type  to  quantify  the 
bead  mass  and  bead  size;  and  (c)  help  parameterize  intra-  and 
interchain  coarse-grained  force-fields  (a  set  of  mathematical 
expressions  describing  various  contributions  to  the  potential 
energy  function  for  the  system  of  interacting  particles,  that  is, 
beads,  in  the  present  case).  A  brief  description  of  the  all-atom 
computational  procedures  employed  in  our  recent  work  [5- 
8, 14]  and  of  the  associated  results  used  in  the  present  work  is 
provided  in  the  following. 

3.1.  PPTA  Coarse-Graining  and  Force-Field  Parameterization. 
Within  the  all-atom  computational  methods,  molecular- 
level  models  of  individual  PPTA  chains  based  on  the  trans 
molecular  conformation  of  the  repeat  unit.  Figure  1(a),  are 
constructed  using  the  Visualizer  [17]  program  from  Accelrys. 
Typically,  PPTA  chains  consisting  of  20  repeat  units  and 
having  the  molecular  weight  of  ca.  4760  Da  are  used. 


of  71,  while,  as  will  be  shown  later,  is  assigned  a  relatively 
large  value. 


Intra-Sheet  Hydrogen-Bond  Term.  Following  our  prior  work 
on  polyurea  [18-20],  the  potential  energy  associated  with 
intra- sheet  interchain  hydrogen  bonding  is  expressed  using 
the  18-9  Lennard-Jones  type  of  function  as 


qihydrogen  _  ^hydrogen 


hydrogen 


18 


-2 


hydrogen 


(3) 


where  is  the  interaction  energy  minimum  and  rep¬ 

resents  the  strength  of  hydrogen-bond  interaction  between 
two  beads,  while  is  the  associated  hydrogen-bonded 

interbead  equilibrium  distance. 


3.1.1.  PPTA  Coarse  Graining.  Since  PPTA  can  be  considered 
as  a  homopolymeric  material  and  in  order  to  be  able  to  handle 
computational  domains  of  a  larger  size,  only  one  type  of  bead 
is  introduced  during  the  coarse-graining  procedure  of  the 
PPTA  chains.  Consequently,  the  bead  mass  was  set  equal  to 
that  of  the  PPTA  repeat  unit  while  the  volume  allotted  to  each 
bead  was  readily  derived  from  the  PPTA- fiber  mass  density. 

3.1.2.  Coarse-Grained  Force-Field  Functions.  The  coarse¬ 
grained  force  field  used  in  the  present  work  includes  four 
(functionally  distinct)  contributions  to  the  system  potential 
energy  (a)  intrachain  bond-stretch  term;  (b)  interchain  two- 
bond-angle  term;  (c)  intra-sheet  hydrogen-bond  term;  and 
(d)  intersheet  van  der  Waals  nonbond  interaction  term.  The 
functional  forms  for  the  four  potential-energy  terms  are  given 
as  follows. 

IntraChain  Bond-Stretch  Term.  Covalent  bonding  within  the 
bead  type  PPTA  chains  is  represented  using  the  following 
form  of  a  harmonic-type  energy  function: 


pStretch  _  ^  ^stretch/  _  stretch 

~  2  y~^o  )  > 


(1) 


where  denotes  the  bond-stiffness  parameter,  r  the 

connector  (i.e.,  bead  bond)  length,  and  the  subscript  o  is  used 
to  denote  the  equilibrium  value  of  the  associated  quantity. 
Due  to  the  aforementioned  high  stretching  stiffness  of  the 
PPTA  molecules,  is  assigned  a  relatively  large  value, 

as  will  be  shown  later. 

IntraChain  Two-Bond-Angle  Term.  An  analogous  harmonic 
function  is  used  to  represent  the  potential  energy  contribu¬ 
tion  associated  with  the  intrachain  two-bond-angle  term  as 


^angle  _  £j^angle^^  _  ^ai 


(2) 


where  denotes  the  bond-stiffness  parameter  and  6  the 
intrachain  interbond  angle.  Due  to  the  aforementioned  high 
bending  rigidity  of  the  PPTA  molecules,  is  set  to  a  value 


Inter  Sheet  van  der  Waals  Nonbond  Interaction  Term.  As  far 
as  the  non-bond  van  der  Waals  interaction  energy  involving 
beads  residing  on  different  PPTA  sheets  is  concerned,  it  is 
represented  using  the  12-6  Lennard-Jones  type  of  function 


as: 


^non-bond  _  ^non-bond 


,,non-bond  \ 


-2 


,^non-bond  \  ^ 


(4) 

where  is  the  interaction  energy  minimum  and 


represents  the  strength  of  nonbond  interaction  between  two 
beads,  while  r^°-bond  ^ 
equilibrium  distance. 


beads,  while  is  the  associated  nonbonded  interbead 


3.1.3.  Coarse-Grained  Force-Field  Parameterization.  Proce¬ 
dures  used  for  parameterization  of  each  of  the  four  parts  of 
the  coarse-grained  force  field  are  presented  in  this  section. 

IntraChain  Bond-Stretch  Term.  Since  the  coarse-grained  for¬ 
mulation  used  in  the  present  work  defines  only  one  type  of 
bead,  the  bond  stretch  force-field  term  contains  only  two 
parameters:  and  ^stretch  determined  by 

setting  it  equal  to  the  equilibrium  bead/bead  separation, 
as  computed  using  all-atom  molecular  statics  calculations 
of  the  minimum  potential-energy  configuration  of  a  single 
infinitely  long,  perfectly  straight  PPTA  chain. 

To  determine  an  equilibrium  molecular-dynam¬ 

ics-based  free-energy  calculation  procedure  is  combined  with 
the  so-called  ''Adaptive  Biasing  Force''  (ABF)  method  [19]. 
During  the  application  of  the  ABF  method,  the  distance 
between  the  centers  of  mass  of  the  (all-atom  represented) 
beads  is  varied  in  small  increments  from  its  initial  (mini¬ 
mum  energy)  value  by  appropriately  changing  the  relative 
position  of  all  constituent  atoms.  This  was  followed  by 
prolonged  (40  ns  simulation  time,  at  a  temperature  of  300  K) 
all-atom  equilibrium  molecular  dynamics  simulations.  The 
trajectory  results  obtained,  averaged  over  the  simulation 
time,  are  used  to  compute  the  mean  force  acting  on  the 
beads,  quantifying  the  tendency  of  the  beads’  connector 
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to  recover  its  equilibrium  length.  To  cancel  this  force,  a 
biasing  force  of  equal  magnitude  and  opposite  sign  is  equally 
partitioned  and  applied  to  all  the  constituent  atoms  in  the 
two  beads.  This  leaves  the  system  under  zero  mean  force  and 
only  the  fluctuating  part  of  the  force  remains,  enabling  the 
computation  of  the  ''equilibrium'  thermodynamic  properties 
of  the  system,  for  example,  mean  potential  energy  at  a  given 
value  of  the  bead  connector  length.  When  the  resulting 
potential  energy  versus  connector  length  results  are  fitted 
using  linear  regression  analysis  in  conjunction  with  (1), 
the  sought-after  stretching  stiffness  parameter  is  obtained. 
The  resulting  functional  relationship  for  bond-stretch  energy 
versus  bonded-bead  separation  is  depicted  in  Figure  5(a). 

IntraChain  Two-Bond-Angle  Term.  The  two-bond-angle 
force-field  term  also  contains  only  two  parameters: 
and  As  mentioned  earlier,  is  set  to  a  value  of  n. 
As  far  as  the  bending  stiffness  parameter  concerned, 

it  is  determined  using  an  ABF-based-procedure  analogous 
to  that  described  for  the  case  of  The  main  difference 

is  that  instead  of  a  perfectly  straight  PPTA  molecule,  a 
molecule  consisting  of  two  perfectly  straight  halves  with 
an  included  angle  6  (of  varying  magnitude)  is  utilized.  The 
resulting  functional  relationship  for  bond-angle  energy 
versus  two -bond- angle  is  depicted  in  Figure  5(b). 

Intra-Sheet  Hydrogen-Bond  Term.  Due  to  the  presence  of 
a  single  bead  type,  the  hydrogen-bond  part  of  the  coarse¬ 
grained  force  field  is  also  characterized  by  two  parameters: 
^hydrogen  ^hydrogen  parameters  Were  determined 

using  the  same  procedure  as  the  ones  described  above,  except 
that  instead  of  a  single  infinitely  long,  perfectly-straight 
chain,  a  wall  of  parallel,  equally  spaced,  infinitely-long, 
perfectly  straight  chains  is  used,  and  the  interchain  spacing 
is  varied.  The  resulting  functional  relationship  for  hydrogen- 
bond  energy  versus  intra-sheet  interchain  separation  distance 
is  depicted  in  Figure  5(c). 

InterSheet  van  der  Waals  Non-Bond  Interaction  Term.  The 
two  non-bond  interaction  parameters  and 

are  determined  using  the  same  procedure,  except  that  an 
infinitely  wide/deep/tall  stack  of  PPTA  sheets  (with  the 
ABABAB. . .  stacking),  each  consisting  of  an  array  of  parallel, 
equally  spaced,  infinitely  long,  perfectly  straight  chains,  is 
used.  The  resulting  functional  relationship  for  nonbond 
energy  versus  intersheet  separation  distance  is  depicted  in 
Figure  5(d). 

Parameterization  Summary.  A  summary  of  the  eight  coarse¬ 
grained  PPTA  force-field  parameters  is  provided  in  Table  1. 


3.2.  Coarse-Grained  Molecular  Method.  All-atom  and  coarse¬ 
grained  molecular  methods  and  tools  are  fundamentally 
identical,  since  both  treat  materials  as  assemblies  of  interact¬ 
ing  particles  and  use  force  fields  to  quantify  various  bonding 
and  nonbonding  interactions.  The  two  main  differences 
between  the  two  methods  are  (a)  the  nature  and  size/mass 
of  the  constituent  interacting  (atoms  versus  beads)  particles; 


Table  1:  Coarse-grained  PPTA  force-field  parameterization  used  in 
the  present  work.  Please  see  (l)-(4)  for  explanation  of  the  symbols. 


Parameter 

Unit 

Value 

j^stretch 

kj/mol-nm^ 

6831.1 

stretch 

nm 

1.31 

j^angle 

kj/mol-rad^ 

9,752 

jangle 

rad 

TC 

^hydrogen 

kj/mol 

22.3 

hydrogen 
'  0 

nm 

1.39 

^non-bond 

kJ/mol 

11.1 

non-bond 

^0 

nm 

1.52 

and  (b)  bead/bead  interactions  are  generally  quite  more 
compliant  (softer)  than  their  atom/atom  counterparts. 

In  general,  a  complete  formulation  of  a  coarse-grained 
molecular  simulation  problem  includes,  at  a  minimum,  the 
following  four  items:  (a)  a  computational  model  consisting 
of  bead-type  chains  and  larger  structures  (e.g.,  sheets,  fibrils, 
fibers,  etc.);  (b)  a  set  of  interbead  force-field  potentials 
which  describe  accurately  various  bonding  and  nonbonding 
bead/bead  interactions;  (c)  a  computational  method(s)  to  be 
used  in  the  simulations;  and  (d)  formulation  of  the  problem 
to  be  analyzed.  These  four  items  are  discussed  in  greater  detail 
in  the  remainder  of  this  section. 

3.2.1.  Computational  Model(s).  While  the  coarse-grained 
computational  approach  allows  a  substantial  increase  in  the 
dimensions  of  the  computational  cell  analyzed,  a  coarse¬ 
grained  analysis  of  a  segment  of  the  PPTA  fiber  (typical 
diameter  on  the  order  of  10  ^m)  would  still  require  the  use  of  a 
prohibitively  large  number  of  beads.  Considering  the  fact  that 
the  majority  of  the  crystalline/topological  defects  are  located 
in  the  core  of  the  PPTA  fibers  and  that  these  defects  play 
an  important  role  in  the  kink-band  nucleation  process,  the 
computational  domain  used  focuses  on  the  core  section  of  the 
PPTA  fibers.  Furthermore,  since  the  transverse  cross-section 
of  the  computational  domain  is  still  small  in  comparison  to 
that  of  the  PPTA-fiber  core,  the  computational  domain  is 
assumed  to  extend  indefinitely  in  the  lateral  directions  (as 
well  as  in  the  axial  direction).  The  rectangular-parallelepiped¬ 
shaped  computational  cell  used  in  the  present  work  is  shown 
in  Figure  6(a).  The  long  edge,  c,  of  the  computational  cell 
is  aligned  with  the  fiber  axis  (which  is,  in  turn,  aligned 
with  the  ^-coordinate  axis).  The  remaining  two  transverse 
edges  (a  and  b)  of  the  cell  are  aligned,  respectively,  with  the 
coordinate  axes  x  and  y.  Typically,  computational  cells  of  the 
size  a  =  24.3  nm,  b  -  24.3  nm,  and  c  =  28.6  nm  are  used.  The 
infinite  extension  of  the  computational  domain  in  the  three 
orthogonal  directions  is  achieved  by  applying  the  periodic 
boundary  conditions  across  the  computational  cell  faces. 

Examination  of  the  transverse  cross-section  of  the  com¬ 
putational  cell  shows  that  the  cell  consists  of  nine  coaxial  fib¬ 
rils  (identified  by  different  colors).  While  each  fibril  contains 
the  same  PPTA  crystal  structure,  the  transverse  orientation 
of  this  structure  differs  among  the  fibrils.  However,  in  all 
fibrils,  the  orientation  of  the  crystal  lattice  is  such  that  the 
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(c)  (d) 

Figure  5:  Coarse-grained  PPTA  force-field  functions  used  in  the  present  work:  (a)  bond-stretch  term;  (b)  two-bond-angle  term;  (c)  hydrogen- 
bond  term;  and  (d)  non-bond  term. 


Figure  6:  Coarse-grained  PPTA  fiber  computational  cell  consisting  of  nine  coaxial  fibrils  (indicated  using  different  colors):  (a)  zoomed-out 
view;  and  (b)  a  close-up  view. 
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sheet  normal  is  orthogonal  to  the  fiber  axis.  The  mean  mass- 
density  within  each  fibril  is  initially  set  to  be  the  same 
(by  properly  adjusting  interchain  distances  within  the  sheets 
and  intersheet  distances)  for  each  fibril.  Before  conducting 
any  analysis/testing  of  the  computational  cell,  the  potential 
energy  of  the  computational  cell  is  minimized  with  respect  to 
the  spatial  coordinates  of  the  constituent  beads.  It  should  be 
noted  that  the  atomic  configuration  displayed  in  Figure  6(a) 
is  the  initial  (before  energy  minimization)  configuration.  A 
close  up  of  this  configuration  is  displayed  in  Figure  6(b).  This 
figure  clearly  shows  the  formation  of  PPTA  sheets  which 
contain  covalent  intrachain  bonds  and  interchain,  intra¬ 
sheet  hydrogen  bonds  (both  depicted  using  the  ball-and-stick 
representation). 

3.2.2.  Force  Fields.  The  force-field  functions  used,  the  proce¬ 
dure  employed  for  their  parameterization  and  the  resulting 
parameters  identification  were  all  covered  in  Section  3.1  and, 
hence,  will  not  be  repeated  here. 

3.2.3.  Computational  Method(s).  Within  the  coarse-grained 
fiber-scale  computational  investigation  carried  out  in  the 
present  work,  two  types  of  computational  methods,  molecu¬ 
lar  statics  (i.e.,  potential-energy  minimization)  and  (equilib¬ 
rium)  molecular  dynamics  (i.e.,  particle-trajectory  analysis) 
are  used.  These  methods  are  implemented  in  Discover  (a 
general-purpose  atomistic  simulation  program)  from  Accel- 
rys  [21],  the  computational  tool  which  was  employed  in 
the  present  work.  Within  the  molecular  statics  approach, 
the  potential  energy  of  the  subject  discrete-particle  material 
system  is  minimized  with  respect  to  the  position  of  the 
particles/beads.  On  the  other  hand,  within  the  molecular 
dynamics  approach,  gradient  of  the  potential  energy  (as 
defined  by  the  force-field  functions)  with  respect  to  the 
particle  positions  is  first  used  to  generate  forces  acting  on 
the  particles,  and  then,  the  associated  Newtons  equations  of 
motion  (for  all  particles)  are  integrated  numerically  in  order 
to  determine  the  temporal  evolution  of  the  particle  positions 
(i.e.,  particle  trajectories).  In  addition,  within  the  equilibrium 
molecular- dynamics  methods,  the  system  under  consider¬ 
ation  is  coupled  to  an  (external)  environment  (a  constant 
temperature  reservoir,  in  the  present  case)  which  ensures 
that  the  system  remains  in  thermodynamic  equilibrium  (i.e., 
the  system  does  not  experience  fluxes  of  thermodynamic 
quantities).  In  the  present  work,  an  NVT  (where  N  is  the 
(fixed)  number  of  atoms,  y,  the  computational  cell  volume 
(also  fixed,  and  set  to  an  arbitrarily  large  value  to  provide  for 
the  isolated-beads  interaction  environment),  and  T  (=300  K) 
is  the  temperature)  rendition  of  the  molecular  dynamics 
method  is  employed. 

As  will  be  shown  in  the  following,  combined  application 
of  the  molecular  statics  and  equilibrium  molecular  dynamics 
methods  enables  the  investigation  of  (a)  axial-compressive 
behavior  of  the  PPTA  fiber  core;  (b)  formation  of  kink 
bands;  and  (c)  elucidation  of  the  role  different  crystallo¬ 
graphic/topological  defects  play  in  the  process  of  kink-band 
formation.  All  the  calculations  in  the  present  work  were 
carried  out  by  coupling  a  user-defined  coarse-grained  force- 
field  database  with  the  Discover  solver  [21].  The  procedure 


employed  to  simulate  PPTA  fiber-fiber  axial  compression 
is  implemented  and  linked  with  Discover  through  the  use 
of  an  input  file  written  in  the  BTCL  script.  This  enabled 
precise  control  of  all  simulation  tasks  such  as  sequential 
contraction  of  the  computational  cell  in  the  longitudinal 
direction  (during  the  simulation  of  the  axial-compression 
test).  In  a  few  cases,  axial- tension  of  the  fibers  which  were 
first  subjected  to  longitudinal  compression  was  investigated 
and,  in  this  case,  loading  is  also  applied  through  the  use  of  an 
input  file  written  in  the  BTCL  language. 

3.2.4.  Problem  Formulation.  As  mentioned  earlier,  the  prob¬ 
lem  analyzed  in  the  present  work  involves  the  mechan¬ 
ical  response  of  PPTA  fibers  subjected  to  longitudinal- 
compression  (and,  in  a  few  cases,  subsequently  subjected  to 
longitudinal  tension).  Typically,  PPTA  fiber  diameter  is  on 
the  order  of  10  pm  and  thus  10^-10^  times  larger  than  the 
transverse  edge  length  of  the  coarse-grained  computational 
cell,  analyzed  in  the  present  work.  As  established  earlier, 
the  computational  cell  used  represents  the  material  located 
within  the  fiber  core  (the  region  of  the  fibers  which  has  the 
highest  probability  for  the  presence  of  flaws/defects  and  thus 
the  highest  potential  for  kink-band  formation).  In  a  typical 
coarse-grained  axial  loading  analysis  of  the  PPTA  fibers,  both 
the  initial  compression  and,  where  applicable,  subsequent 
longitudinal  tension  were  applied  in  small  increments  and 
each  incremental  deformation  is  followed  by  the  compu¬ 
tational  cell  energy-minimization  and  extended  (ca.  40  ns) 
molecular- dynamics  simulation  runs. 

3.2.5.  Postprocessing  Data-Reduction  Analysis.  The  coarse¬ 
grained  fiber-level  particle- trajectory  computational  results 
obtained  are  mainly  used  to  reveal  the  response  of  sin¬ 
gle  PPTA  fibers  to  uniaxial  compression  (including  the 
microstructural/topological  evolution  processes  such  as  kink 
band  formation  and  fibrillation).  Fiber  axial  compression  is 
simulated  by  subjecting  the  unit  cell  to  longitudinal  compres¬ 
sive  deformation  in  the  ^-direction  while  allowing  for  the 
operation  of  the  Poissons  effects  in  the  x-  and  y-directions.  In 
this  case,  except  for  graphical  display,  no  data  postprocessing 
procedure  had  to  be  employed.  However,  in  a  few  simulations, 
the  effect  of  these  processes  on  the  residual  axial-tensile 
strength  of  the  fiber  was  investigated.  In  these  cases,  (a) 
following  uniaxial  compression  and  the  formation  of  kink 
bands/fibrillation,  the  unit  cell  is  subjected  to  longitudinal 
tensile  deformation  while  allowing  for  the  operation  of  the 
Poissons  effects  in  the  transverse  directions;  and  (b)  the  tra¬ 
jectory  data  had  to  be  postprocessed  to  extract  and  quantify 
the  effect  sought  after.  During  the  application  of  the  tensile 
deformation,  the  procedure  described  in  our  recent  work  [18] 
was  utilized  to  compute  the  normal  axial  stress  and  its  rate  of 
change  increase  with  an  increase  in  the  axial  tensile  strain. 
The  longitudinal-tensile  strength  of  the  single  PPTA  fibril 
is  then  set  equal  to  the  value  of  the  corresponding  normal 
stress  at  which  the  rate  of  stress  increase  with  increase  in 
strain  begins  to  decline  ''appreciably'  (denoting  the  onset  of 
inelastic- deformation  and/or  damage -initiation  processes). 
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4.  Results  and  Discussion 

In  this  section,  the  results  of  the  coarse-grained  fiber-level 
computational  analyses  are  used  to  assess  (a)  axial  compres¬ 
sive  response  of  single  PPTA  fibers;  (b)  kink-band  formation 
and  fibrillation  processes;  (c)  the  effect  of  prior  longitudinal 
compression  on  the  residual  axial  tensile  strength  of  single 
PPTA  fibers;  and  (d)  the  role  of  various  micro  structural  and 
topological  defects  on  the  kink-band  formation/fibrillation 
propensity  of  the  PPTA  fibers  and  the  effect  of  these  localized 
deformation/damage  processes  on  the  residual  longitudinal- 
tensile  strength  of  the  fiber. 

Since  the  work  presented  in  this  section  considers  PPTA 
fibers  which  are  initially  either  perfectly  crystalline  or  crys¬ 
talline  but  contain  various  defects,  fibers  which  are  subjected 
to  tensile  loading  directly  or  after  being  subjected  to  a 
preludial  axial  compression  and  attention  is  paid  to  both 
defect- induced  and  compression- induced  changes  in  the 
material  microstructure  and  the  longitudinal  tensile  strength, 
the  corresponding  results  are  presented  and  discussed  within 
separate  sections  (each  dealing  with  a  particular  case  of  the 
fiber  crystallinity  and  prior  axial  compression). 

4.1.  Longitudinal-Tensile  Strength  of  Defect-Free  Virgin  PPTA 

Fibers.  The  procedure  for  calculation  of  the  fiber  longitudi¬ 
nal-tensile  strength  was  described  in  Section  3.2.  This  proce¬ 
dure,  in  the  case  of  defect-free  virgin  PPTA,  yielded  a  fairly 
repeatable  value  of  the  fibril  longitudinal-tensile  strength 
of  approximately  6.3  GPa.  Close  examination  of  the  fiber 
microstructure  at  the  onset  of  inelastic  deformation/damage 
initiation  revealed  that,  in  this  case,  the  longitudinal-tensile 
strength  is  controlled  by  the  onset  of  pronounced  degrada¬ 
tion  of  interfibril  bonding  as  well  as  interchain,  intra-sheet 
bonding.  It  should  be  noted  that  the  term  is  used 

to  denote  the  (either  perfectly  crystalline  or  defective)  PPTA 
fibers  which  were  not  subjected  to  axial  compression  prior  to 
longitudinal  tension. 

4.2.  Microstructure  Evolution  during  Axial  Compression  of 
Defect-Free  PPTA  Fibers.  Evolution  of  the  coarse-grained 
microstructure,  during  axial  compression,  within  the  com¬ 
putational  cell  containing  (initially)  defect-free  virgin  PPTA 
fibers,  is  depicted  in  Eigures  7(a) -7(d).  The  longitudinal 
stretch  values  associated  with  Eigures  7(a)-7(d)  are  (a)  0.975; 
(b)  0.950;  (c)  0.925;  and  (d)  0.900.  The  lateral  stretches, 
in  each  case,  take  on  the  values  required  by  the  zero- 
lateral-stress  condition.  Examination  of  the  results  displayed 
in  these  figures  clearly  reveals  formation  and  subsequent 
growth/thickening  of  two  counter-balancing  nonparallel 
shear  bands. 

The  shear  bands  formed  are  found  to  be  of  a  reversible 
nature,  that  is,  they  vanish  upon  removal  of  the  axial  com¬ 
pressive  loading.  Detailed  examination  of  the  microstructure 
obtained  after  removal  of  the  axial-compressive  loading 
revealed  that  this  microstructure  is  very  similar  to  that 
observed  in  the  defect-free  virgin  PPTA  fiber.  In  other  words, 
preludial  axial  compression  in  the  case  of  defect-free  PPTA 
fibers  does  not  cause  any  significant  permanent  change  in  the 
fiber  microstructure. 


4.3.  Longitudinal-Tensile  Strength  of  Defect-Free  Precom¬ 
pressed  PPTA  Fibers.  Application  of  the  aforementioned 
longitudinal-tensile-strength  determination  procedure  to  the 
defect-free  precompressed  PPTA  fibers  revealed  that  this 
strength  is  practically  identical  to  its  counterpart  (6.3  GPa) 
obtained  in  the  case  of  defect-free  virgin  PPTA  fibers.  This 
finding  is  fully  consistent  with  the  aforementioned  obser¬ 
vation  that  the  preludial  compression  (at  least  up  to  the 
compressive  strain  of  ln(0.9)  =  -0.105)  does  not  introduce 
any  significant  permanent  changes  in  the  fiber  material 
microstructure. 

4.4.  Microstructure  of  Defective  Virgin  PPTA  Fibers.  In  our 
recent  work  [5-8,  14],  a  detailed  overview  was  provided 
with  respect  to  the  character  of  the  main  microstruc- 
tural/ topological  defects  found  in  PPTA  fibers,  their  clas¬ 
sification,  their  origin  in  relation  to  the  PPTA  synthe¬ 
sis/processing  conditions,  and  their  (estimated)  typical  con¬ 
centrations.  Consequently,  a  similar  detailed  overview  will 
not  be  provided  here.  However,  it  should  be  established 
that  among  the  isolated  and  clustered  point-type  defects, 
the  following  were  identified  as  having  the  largest  effect  on 
the  PPTA- fiber  properties  (in  particular,  on  the  longitudinal- 
tensile  strength):  (a)  chain-ends;  (b)  side  groups;  and  (c)  voids 
and  interstitials.  In  addition,  planar-type  defects,  commonly 
referred  to  as  defect  bands,  were  also  found  to  have  a  major 
effect  on  the  fiber  strength.  In  the  present  work,  only  the 
effect  of  (isolated  or  clustered)  chain  ends  or  side  groups  was 
investigated.  Chain  ends  are  the  defects  which  stem  from  the 
fact  that  some  of  the  fibril  chains  do  not  run  continuously 
along  the  fiber  length  but  rather  terminate  (at  one  or  both 
ends)  within  the  fiber.  As  far  as  the  side-group  defects  are 
concerned,  they  are  associated  with  the  presence  of  radicals, 
for  example,  -SO3H,  -S03~Na'^,  between  the  PPTA  sheets. 
As  shown  in  our  previous  work  [5-8, 14] ,  the  presence  of  these 
two  types  of  defects  and  their  clusters  within  the  (otherwise) 
perfect  PPTA  fibers  causes  local  changes  in  the  material 
microstructure.  These  microstructural  changes,  as  will  be 
shown  in  the  following,  affect  the  fiber  longitudinal-tensile 
strength. 

4.5.  Longitudinal-Tensile  Strength  of  Defective  Virgin  PPTA 
Fibers.  Application  of  the  longitudinal-tensile-strength  de¬ 
termination  procedure  clearly  revealed  that  the  fiber  strength 
is  highly  sensitive  to  the  presence,  type,  and  cluster  size  of  the 
two  types  of  defects  analyzed.  Due  to  the  stochastic  nature 
of  the  spatial  distribution  of  the  defect  concentration  within 
the  fiber,  the  longitudinal-tensile  strength  of  the  fibers  has 
also  been  found  to  be  a  stochastic  quantity.  This  finding  is  in 
full  agreement  with  several  experimental  reports  (e.g.,  [22- 
24]).  In  our  recent  work  [6],  a  procedure  was  introduced  for 
the  derivation  and  parameterization  of  a  probability  density 
function  for  the  longitudinal-tensile  strength  of  PPTA  fibers. 
This  procedure  will  be  used  in  the  present  work  (without 
reviewing  the  details).  When  applying  this  procedure,  (iso¬ 
lated  and  clustered)  defect  concentrations  were  set  to  their 
values  associated  with  the  prototypical  PPTA  fiber  fabrication 
conditions  [6]. 
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Figure  7:  Coarse-grained  microstructure,  during  axial  compression,  within  the  computational  cell  containing  (initially)  defect-free  virgin 
PPTA  fiber  at  the  longitudinal  stretch  values  of  (a)  0.975;  (b)  0.950;  (c)  0.925;  and  (d)  0.900. 


Application  of  this  procedure  to  the  results  obtained 
in  the  present  work  yielded  the  virgin-defective-crystal 
PPTA-fiber  longitudinal-tensile-strength  probability  density 
and  (complementary)  cumulative  density  functions  depicted, 
respectively,  in  Figures  8(a)  and  8(b).  It  should  be  noted  that 
the  probability  density  function  in  Figure  8(a)  is  represented 
by  its  corresponding  relative  frequency,  while  the  comple¬ 
mentary  cumulative  density  function  displayed  in  Figure  8(b) 
is  obtained  by  integrating  the  relative  frequency  from  oo  to 
the  given  value  of  the  fiber  axial  strength.  The  mean  value 
of  the  virgin-defective-PPTA  longitudinal  tensile  strength, 
computed  as  the  first  moment  of  the  probability  density 
function  displayed  in  Figure  8(a),  is  found  to  be  ~5.0  GPa. 
This  value  is  in  reasonably  good  agreement  with  the  measured 
PPTA-fiber  longitudinal-tensile  strength  (typically  found  to 
be  in  a  3-5  GPa  range  (e.g.,  [20,  23,  25,  26]).  This  find¬ 
ing  suggests  that  the  present  computational  procedure  can 
reasonably  well  account  for  the  experimental  observations 


regarding  the  PPTA  fiber  longitudinal-tensile  strength.  The 
cumulative  distribution  displayed  in  Figure  8(b)  defines  the 
probability  that  the  PPTA  fiber  longitudinal-tensile  strength 
is  not  lower  than  a  given  (lower  bound)  value,  and  this  type  of 
function  is  often  used  in  probabilistic/reliability-based  design 
analyses  [27].  In  such  analyses,  functional  constraints  cannot 
be  met  with  absolute  certainty  but  rather  with  an  acceptably 
low  level  of  failure. 

4.6.  Microstructure  Evolution  during  Axial  Compression  of 
Defective  PPTA  Fibers.  Introduction  of  chain-end  isolated 
or  clustered  defects  into  the  PPTA  fiber  did  not  change 
the  nature  of  the  microstructure  evolution  during  longitu¬ 
dinal  compression.  In  other  words,  axial  compression  was 
accompanied  by  the  formation  and  growth  of  counterbal¬ 
ancing  kink  bands.  However,  kink  bands  are  found  to  occur 
at  a  lower  level  of  axial- compressive  strain  (revealing  the 
stimulating  effect  that  these  defects  have  on  the  kink-band 
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Figure  8:  (a)  Longitudinal- tensile- strength  relative  frequency  (i.e.,  probability  density  function);  and  (b)  the  corresponding  complementary 
cumulative  density  function  for  defective  virgin  PPTA-fibers. 


formation).  Since  the  fiber  microstructures,  at  different  stages 
of  axial  compression,  obtained  in  this  case  are  very  similar  to 
those  displayed  in  Figures  7(a) -7(d),  they  are  not  shown  here. 
Another  difference  between  the  axial-compressive  behavior 
of  the  defect-free  and  the  defective  PPTA  fibers  is  that,  in 
some  cases,  removal  of  the  axial- compressive  loading  did  not 
result  in  complete  disappearance  of  the  kink  bands. 

When  isolated  or  clustered  side-group  defects  were 
introduced  into  the  PPTA  fiber,  axial  compression  was  not 
accompanied  by  the  formation  of  well-defined  kink  bands. 
Rather,  microstructural  changes  resembling  fibrillations  were 
observed  [28,  29].  An  example  of  the  microstructural  evolu¬ 
tion  during  axial  compression  of  PPTA  fibers  containing  side- 
group  defects  is  depicted  in  Figures  9 (a) -9(d).  The  longitu¬ 
dinal  stretch  values  associated  with  Figures  9(a)-9(d)  are  (a) 
0.975;  (b)  0.950;  (c)  0.925;  and  (d)  0.900.  The  lateral  stretches, 
in  each  case,  take  on  the  values  required  by  the  zero -lateral- 
stress  condition.  Examination  of  the  results  displayed  in  these 
figures  reveals  the  formation  of  bulges  (an  indication  of  the 
operation  of  fibrillation  processes)  at  the  computational- cell 
lateral  faces.  The  height  and  the  extent  of  these  bulges  increase 
in  the  course  of  axial  compression.  In  most  cases  analyzed, 
microstructural  changes  experienced,  during  preludial  com¬ 
pression,  by  the  PPTA  fiber  containing  side-group  defects 
were  permanent.  That  is,  with  the  exception  of  minor  elastic- 
relaxation  processes,  removal  of  the  compressive  loading  did 
not  restore  the  original/uncompressed  microstructure. 

4.7.  Longitudinal-Tensile  Strength  of  Defective  Precompressed 
PPTA  Fibers.  In  the  case  of  PPTA  fibers  containing  chain-end 
defects,  application  of  the  prior  axial  compression  was  found 
not  to  measurably  change  the  longitudinal-tensile-strength 
probability  density  and  cumulative  density  functions  (the 
results  not  shown  for  brevity)  relative  to  their  counterparts 


observed  in  the  uncompressed  PPTA  fibers  containing  the 
same  type/number  of  defects.  This  finding  is  fully  consistent 
with  the  fact  that,  in  most  cases,  kink  bands  formed  during 
axial  compression  of  the  PPTA  fibers  containing  chain-end 
defects  vanish  upon  the  removal  of  the  compressive  loading 
with  no  major  residual  changes  in  the  fiber  microstructure. 

In  sharp  contrast  to  the  finding  just  described,  in  the 
case  of  PPTA  fibers  containing  side-group  defects,  preludial 
axial  compression  is  found  to  change  the  longitudinal-tensile 
strength.  Typical  results  obtained  in  this  case  are  presented 
in  Figures  10(a),  and  10(b).  The  results  displayed  in  Figures 
10(a)  and  10(b)  pertain  to  the  probability  density  and  the 
cumulative  distribution,  respectively.  A  comparison  of  the 
corresponding  results  displayed  in  Figures  8(a),  8(b),  10(a) 
and  10(b)  reveals  that  the  presence  of  side-group  defects 
and  the  application  of  preludial  axial  compression  result 
in  a  reduction  in  the  mean  value  (4.86  GPa)  of  the  fiber 
longitudinal-tensile  strength.  In  addition,  the  distribution 
standard  deviation  is  increased  from  -0.40  GPa  (for  the  case 
of  the  same  fibers  which  were  not  subjected  to  preludial 
axial  compression)  to  -0.52  GPa  (in  the  present  case).  This 
finding  suggests  that,  in  the  case  of  precompressed  fibers, 
there  is  a  higher  probability  for  finding  fibers  with  relatively 
low  longitudinal-tensile  strength. 

The  results  presented  in  this  section  suggest  that  the 
nature  (i.e.,  chain  end  versus  side  group)  of  the  defects, 
as  well  as  the  defect  concentration/cluster- size,  may  play 
a  key  role  in  the  effect  of  preludial  axial  compression  on 
the  residual  longitudinal-tensile  strength  of  the  PPTA  fibers. 
That  is,  preludial  axial  compression  has  substantially  smaller 
detrimental  effect  on  the  fiber  longitudinal-tensile  strength  in 
the  case  of  chain -end  defects.  In  sharp  contrast,  this  negative 
effect  can  become  significant  when  PPTA  fibers  contain  side- 
group  defects. 
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Figure  9:  Coarse-grained  microstructure,  during  axial  compression,  within  the  computational  cell  containing  (initially)  defective  virgin 
PPTA  fiber  at  the  longitudinal  stretch  values  of:  (a)  0.975;  (b)  0.950;  (c)  0.925;  and  (d)  0.900.  The  fiber  within  the  computational  cell  contains 
a  cluster  of  side-group  defects. 


5.  Summary  and  Conclusions 

Based  on  the  results  obtained  in  the  present  work,  the 
following  summary  remarks  and  main  conclusions  can  be 
drawn. 

(1)  The  effect  of  prior  axial  compression  on  the  resid¬ 
ual  longitudinal-tensile  strength  of  defect-free  and 
defective  p-phenylene  terephthalamide  (PPTA)  fibers 
is  investigated  computationally  using  coarse-grained 
molecular  statics  and  dynamics  simulation  methods. 
To  construct  the  coarse-grained  PPTA-fiber  model, 
the  results  of  a  companion  all-atom  molecular-level 
computational  analysis  are  utilized. 

(2)  In  the  case  of  defect-free  PPTA  fibers,  the  longitu¬ 
dinal-tensile  strength  was  found  to  be  essentially  a 
deterministic  quantity.  In  sharp  contrast,  in  defective 
PPTA  fibers,  the  same  quantity  was  found  to  be  of 
a  stochastic  character  and  to  possess  a  characteristic 
probability  density  function. 

(3)  Preludial  compression  was  found  to  alter  the  low- 
order  moments  of  this  distribution  function,  and  the 


extent  of  these  alterations  was  found  to  be  dependent 
on  the  defect  type  and  concentration.  Specifically, 
chain- end  defects  are  found  to  make  this  alteration 
small,  while  side  groups  are  observed  to  increase  the 
extent  of  these  alterations.  In  the  case  of  PPTA  fibers 
containing  side-group  defects,  preludial  compression 
has  been  found  to  reduce  the  mean  value  of  the 
fiber  longitudinal-tensile  strength  and  to  increase  the 
variance/standard-deviation  of  the  associated  prob¬ 
ability  density  function.  The  extent  of  these  changes 
was  observed  to  increase  with  the  extent  of  preludial 
axial  compression. 

(4)  The  results  obtained  clearly  reveal  that  the  quality 
control  related  to  the  precursor  materials,  synthe¬ 
sis  conditions,  and  fiber-spinning  process  must  be 
given  particular  attention,  since  the  accompanying 
inadequacies  may  result  in  PPTA  fibers  with  inferior 
mechanical  properties. 

(5)  The  results  also  point  out  that  various  textile¬ 
manufacturing  processes  used  in  the  fabrication  of 
high-performance  protective  systems  may  degrade 
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Figure  10:  (a)  Longitudinal-tensile- strength  relative  frequency  (i.e.,  probability  density  function)  and  (b)  the  corresponding  complementary 
cumulative  density  function  for  PPTA-fibers  containing  side-group  defects  after  preludial  axial  compression  to  a  maximum  strain  of  -0.105. 


the  properties  of  PPTA  fibers  relative  to  their  as-syn¬ 
thesized/as-spun  levels. 
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